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Abstract 

The framework of emergent gravity arising from Yang-Mills matrix models is 
developed further, for general noncommutative branes embedded in M D . The 
effective metric on the brane turns out to have a universal form reminiscent of 
the open string metric, depending on the dynamical Poisson structure and the 
embedding metric in IR D . A covariant form of the tree- level equations of motion 
is derived, and the Newtonian limit is discussed. This points to the necessity 
of branes in higher dimensions. The quantization is discussed qualitatively, 
which singles out the IKKT model as a prime candidate for a quantum theory 
of gravity coupled to matter. The Planck scale is then identified with the scale 
of N = 4 SUSY breaking. A mechanism for avoiding the cosmological constant 
problem is exhibited. 
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1 Introduction 

The notion of space-time which underlies the presently accepted models of fundamental 
matter and interactions goes back to Einstein. Space-time is modeled by a 4-dimensional 
manifold, whose geometry is determined by a metric with Lorentzian signature. This notion 
escaped the quantum revolution essentially unchanged, even though Quantum Mechanics 
combined with General Relativity strongly suggests a "foam-like" or quantum structure at 
the Planck scale. While some kind of quantum structure of space-time indeed arises e.g. 



in string theory or loop quantum gravity, a satisfactory understanding is still missing. 



A different approach to this problem has been pursued in recent years, starting with 
some explicit quantization of space-time and attempting to construct physical models on 
such a background. The classical space-time M 4 is replaced by a quantized or "noncom- 
mutative" (NC) space, where the coordinate functions x^ satisfy nontrivial commutation 
relations such as [x^, x v \ = i9^ v '. This leads to non-commutative field theory, see e.g. [0-0]. 
At the semi-classical level, these commutation relations determine a Poisson structure 9^ 
on space-time, which is fixed by construction. However, since quantized spaces are expected 
to arise from quantum gravity, it seems more appropriate to consider a dynamical Poisson 
structure at the semi-classical level. A straightforward generalization of General Relativity 
is then inappropriate; indeed any quantum structure of space-time rules out classical intu- 
itive principles. Rather, one should look for simple models of dynamical noncommutative 
(or Poisson) spaces, with the hope that they will effectively incorporate gravity. 

Such models are indeed available and known as Matrix-models of Yang-Mills type. 
They have the form S = Tr[X a ,X b ][X a , X b ]S aa iSw + where indices run from 1 to 
D. It is well known that these models admit noncommutative spaces ("NC branes") as 
solutions, such as the Moyal-Weyl quantum plane 1^; see e.g. |T3Hl9|. However, most 



of the work up to now is focused on special NC branes with a high degree of symmetry. 
For generic NC spaces with non-constant 9 ,xv {x), it was shown in Q that the kinetic term 
for any "field" coupled to the D = 4 matrix model is governed by an effective metric 
G ab (x) = p9 aa (x)9 bb {x)8 a >y, including nonabelian gauge fields. This nicely explains the 
observed relation in || between NC U(l) gauge fields and gravitational degrees of freedom, 
see also |6|-[8| for related work. Since this effective metric is dynamical, these YM Matrix 
Models contain effectively some version of gravity, thus realizing the idea that gravity 
should emerge from NC gauge theory As argued in [Q, an effective action for gravity 

is induced upon quantization, with the remarkable feature that the "would-be cosmological 
term" decouples from the model due to the constrained class of metrics. This makes the 
mechanism of induced gravity feasible at the quantum level, and suggests that the Newton 
constant resp. the Planck scale is related to an effective UV-cutoff of the model. A detailed 



analysis taking into account UV/IR mixing |10j and fermions ]TT) singles out the N = A 



supersymmetric extensions of the model, where such a cutoff is given by the scale of iV = 4 



SUSY breaking. This amounts to D = 10, which is nothing but the IKKT model [ p"2|j , 
originally proposed as a nonperturbative definition of IIB string theoryf]. 

In the present paper, we develop the framework for emergent gravity on general NC 
branes with nontrivial embedding in M D . This works out very naturally, leading to a 
simple generalization of the effective metric which is strongly reminiscent of the open string 



1 As such, the presence of gravity in this model is expected and to some extent verified, cf. p~3|- 
|22|,|4],|5|. However, what is usually considered are effects of D=10 (super)gravity, modeled by interactions 
of separated "D-objects", represented by block-matrices. In contrast, emergent NC gravity describes 
interactions within (generic) NC branes in this model. Evidence for gravity on simple NC branes was 
obtained previously in p4,pr|. 
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metric [p7 , involving the general Poisson tensor and the embedding metric. We establish in 



Section |]the relevant geometry, find the semi-classical form of the bare matrix- model action 
for general NC branes in W D , and obtain covariant equations of motion. This generalizes 
the well-known case of flat or highly symmetric branes to the generic case, and shows how 
the would-be U(l) gauge field is absorbed in the effective metric on the brane. In Section 
|3], the Newtonian limit of emergent gravity is studied in detail. It turns out that even 
though it is possible to reproduce the Newtonian potential for general mass distributions, 
the relativistic corrections are in general not correctly reproduced in D = 4 matrix models. 
This provides one motivation to consider general branes embedded in higher- dimensional 
matrix models, which admit a much richer class of geometries and promise to overcome this 
problem. The compactification of higher-dimensional NC branes is described in Section E]4 
with the example of fuzzy spheres in extra dimensions. 

Higher dimensions, more precisely D — 10 resp. N = 4 SUSY also appears to be 
required by consistency at the quantum level. From the point of view of emergent gravity, 
this condition arises as a result of UV/IR mixing in NC gauge theory. This is discussed 
qualitatively in Section ^]6| along the lines of [|J , leading to an induced gravity action. In 



Section |27| some differences to General Relativity are discussed, most notably the presence 
of intrinsic scales and preferred coordinates, as well as the different role of the "would-be 
cosmological constant term" . As an illustration of the formalism, we also give a (unphysical) 
solution of the bare equations of motion in Section Finally, a matrix version of a conserved 
energy-momentum tensor is derived. 

The results of this paper provide a rich framework for the search of realistic solutions of 
emergent NC gravity. The main missing piece is the analog of the Schwarzschild solution, 
which is nontrivial because the quantum effective action at least at one loop must be taken 
into account. But in any case, it is clear that these models do contain a version of gravity 
in an intrinsically noncommutative way, and they have a good chance to be well-defined at 
the quantum level at least for the IKKT model. This certainly provides motivation for a 
thorough investigation. 



2 The Matrix Model 

Consider the matrix model with action 

S YM = -Tr[X^X v ][X»\X v ']g^g wli (1) 

for 

fiV = <W or SV = Vw? ( 2 ) 
in the Euclidean resp. Minkowski case. The "covariant coordinates" X^, fx = 1,2,3,4 

are hermitian matrices or operators acting on some Hilbert space "H. We will denote the 
commutator of 2 matrices as 

[X^X v ] = i6^ (3) 
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so that 9^ v G L(9C) are antihermitianQ matrices, which are not assumed to be proportional 
to 1%. We focus here on configurations X M which can be interpreted as quantizations of 
coordinate functions x M on a Poisson manifold (Ai, 9 fiU (x)) with general Poisson structure 
9 fJ,u (x). This defines the geometrical background under consideration, and conversely es- 
sentially any Poisson manifold provides (locally) a possible background X^ More 
formally, this means that there is an isomorphism of vector spaces 

C(M) ->• A C L(H) 

f(x) ^ f(X) (4) 

i{f,g] ^ [f,g} + 0(9 2 ) 

Here C(M.) denotes some space of functions on .M, and A is interpreted as quantized 
algebra of functions^ on At. This allows to replace [f(X),g(X)] — > i{f(x),g(x)} to leading 
order in 9^. In particular, we can then write 

[X»,f{X)]~W"( x )±.f( x ) (5) 

which will be used throughout this paper, denoting with ~ the leading contribution in a 
semi-classical expansion in powers of 9^ u . 

In order to derive the effective metric on Ai, let us now consider a scalar field coupled 
to the matrix model (|I|). The only possibility to write down kinetic terms for matter fields 
is through commutators [X^, $] using (|5[). Thus consider the action S = Sym + where 

S[«] = -Trg m ,[X»,$\\X l ''A 

~ jkJ ^ XP(X)G " (X) ^ (X) ^ (X) - <6) 

Here 

cr{x) = e^{x)r / {x)g^ v , (7) 

is interpreted as metric on Ai in x coordinates. We will assume in this paper that 9 fJ, "(x) is 
nondegenerate. Then the symplectic measure on (Ai, 9 fil/ (x)) is given by the scalar density 

= I^WI 172 = IG^WW = AW > (8) 

which can be interpreted as "local" non-commutative scale Atvc- m the preferred x coor- 
dinates characterized by (fj), p(x) coincides with the dimensionless scalar function 



\G» v {x)\^ _ \9-l{x)\^ 



(9) 



2 in contrast to the conventions in [Q 

3 Roughly speaking A is the algebra generated by X^, but technically one usually considers some sub- 
algebra corresponding to well-behaved functions. 
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The action can now be written in a covariant manner as 



Sm = ?2^a / d 4 xG^(x) d^(x)dMx) = tA^ J d 4 x $(x)A e $(x) , (10) 

Here is the Laplacian for the metric [|j] 

(x) = | | 1/4 (x) = e- CT (x) , 
\& v \ = 1, (11) 

which is unimodular in the preferred x M coordinates. By definition, G^ v (x) is the effective 
metric for the scalar field. Because it enters in the kinetic term for any matter coupled to 
the matrix model, it plays the role of a gravitational metric 

ds 2 = G^(x) dx^dx v . (12) 

Up to certain density factors, this also applies to nonabelian gauge fields as shown in [Q and 
for fermions [[II] . Therefore the Poisson manifold under consideration naturally acquires 
a metric structure (Ai, 6^ v {x), (5 MI/ (x)), which is determined by the Poisson structure and 
the flat background metric g^ u . 



Equations of motion. The basic matrix model action (|I|) leads to the e.o.m. for 

p,p',r']] ffwi , = o. (13) 

This can be written in the semi-classical limit as O^d^^^g^ = 0, or 

G»(x)a i 0^ = O. (14) 

These equations are not covariant, they are valid only in the coordinates x M where the 
"background metric" in the matrix model is either 5^ or rj^. As shown in Appendix 
B, these equations of motion can be written in a covariant manner as 



G^(x)V 7 (e CT =e-°G„ v e^d 1 r } {x 



where 



fl/l' VV 



(15) 



(16) 



and V denotes the Levi-Civita connection with respect to the metric G^ . Note that the 
"background" metric g^ u is absorbed completely. (|1^) can be written as 



^(x)V^- 1 



-G 7r? (x) 9-}d j( j + e- 2a G^ d^x) 



(17) 



which has the form of covariant Maxwell equations with source. The obvious advantage 
of this covariant form of the equations of motion is that we can now use any adapted 
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coordinates, in particular rotation-symmetric ones etc. This should help to find solutions. 
Nevertheless, this should not obscure the fact that the underlying matrix model is not 
invariant under diffeomorphisms: the background metric g^ v is constant, and there is no 
obvious way to transform it at the level of the matrix model. Only in the semi-classical 
limit we can allow general coordinates and rewrite things in a coordinate independent way, 
at the expense of introducing a flat background metric g^ v . 

In principle of course, the equation of motion for X^ is modified due to the presence 
of the scalar field. However for small coupling or energy, we can presumably neglect this 
back-reaction of matter on the geometry. It will be taken into account in section ( |2.1| ). 

The equation of motion for the scalar field <fi are 

= [X", \X\ <j>]]g^ ~ e^'d^'d u 4)9,u ■ (18) 

As shown in Appendix A, this can be written as 

= {GTdpdv - V"dJ<t> = (19) 

where = G vri Y^ and T% are the Christoffel symbols of G^. This follows also immedi- 



ately from the covariant form (10) of the scalar action. We will show moreover in Appendix 
A that in the preferred x M coordinates defined by the matrix model, the equation of motion 
(14]) for X^ resp. 9~l is equivalent to the non-covariant equation ( |181|) 



f M = . (20) 

In these coordinates, the equation of motion for <fi takes the simple form G^d^d^ = 0, for 
on-shell geometries. 

The semi-classical form of the matrix model action ([I]) is 

Sym = ttA^- / d 4 x p(x)r](x) . (21) 



(2k) 



The equation of motion (|15D can be derived directly from this action. We will give this 
derivation in the next section, in the context of general branes embedded in IR D . 



2.1 Noncommutative branes and extra dimensions 

Let us discuss scalar matter from the point of view of extra dimensions. Recall that e.g. 
the action for a scalar field is given by additional terms of the type 

Tr^^X",^. (22) 

The combined action can be interpreted as matrix model with extra dimensions, where one 
coordinate denoted as is a function of the other 4 coordinates. Therefore we consider 
more generally 

S Y M = -Tr[X a ,X b ][X a ',X b ']r ]aa ,r ]w , (23) 
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for hermitian matrices or operators X a } a = 1, D acting on some Hilbert space "H. To 
avoid a proliferation of symbols we fix the background to have the Minkowski metric; the 
Euclidean case is completely parallel, replacing r] a f, with 5 a b. A scalar field can therefore 
be interpreted as defining an embedding of a 4-dimensional manifold (a "3-brane") in a 
higher-dimensional space. This naturally suggests to consider a higher- dimensional version 
of the Yang-Mills matrix model, such as the IKKT model in 10 dimensions. 

We want to consider general 2n -dimensional noncommutative spaces M. 2 ) 11 C M. D (a 
2n — 1 brane) in D dimensions. We correspondingly split the matrices as 

X a = // = l,...,2n, i = l,...,D-2n. (24) 

The basic example is a flat embedding of a 4-dimensional NC background with 

[X»,X V ] = i9»\ /x,z/=l,...,4, 

<f? = 0, z = l, ...,£> -4 (25) 

where X^ generates a 4-dimensional NC brane Aig. Then the discussion of the previous 
section applies, and fluctuations of <p l (x) can be interpreted as scalar fields on M.\. More 
generally, we can interpret <p l (x) as defining the embedding of a 2n - dimensional submani- 
fold Ai 2n C M D , equipped with a nontrivial induced metric. The support ( "D-dimensional 
spectrum") of X a ~ x a will then be concentrated on Ai 2n C ~R D in the semi-classical limit. 
Expressing the 4> % in terms of X M , we obtain 

[cj>\ ~ iO^drfdJ = ie^(f) (26) 

in the semi-classical limit. This involves only the components fj, — 1, 2n of the antisym- 
metric tensor [X a ,X fe ] ~ i$ ab (x), which has rank 2n in this case. Here 

e » : = -ifX**,.] ~ O^du (27) 

are derivations, which span the tangent space of M. 2n C M. D . They will define a preferred 
frame below. We can then interpret 

[X»,X v }~i9^(x) (28) 

as Poisson structure on Ai 2n (assuming that it is no-degenerate), noting that the Jacobi 
identity is trivially satisfied. This is the Poisson structure on Ai 2n whose quantization 
is given by the matrices X M , jj, — 1, ...,2n, interpreted as quantization of the coordinate 
functions x^ 1 on M 2n . Conversely, any 6^ u {x) (|28|) can be (locally) quantized, and provides 
together with arbitrary (j) l (x) a quantization of Ai 2n C M 15 as described above. Note that 
this Poisson structure is defined intrinsically by the configurations of the matrix model, 
independent of the choice^ made in fl24]) . Assuming that 6^ v (x) is non-degenerate, we 

4 For generic embeddings, the separation (|24| ) is arbitrary, and we are free to choose different 2n compo- 
nents among the {X a } as generators of tangential vector fields. This is a particular change of coordinates 
on A4 2n , which from the field theory point of view corresponds to a remarkable transformation exchang- 
ing fields with coordinates, reminiscent of T-duality in string theory. In any case, note that iQ^ v is not 
naturally a pull-back of some non-degenerate Poisson or symplectic structure on MP . 
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denote its inverse matrix with 

0~l{x) , (29) 

which defines a symplectic form on Ai 2n . Finally, the trace is again given semi-classically 
by the volume of the symplectic form, 



(2n) n Trf~ J d' n xp(x)f (30) 

with p = (det^J) 1 / 2 generalizing (|[). 

We are now in a position to extract the semi-classical limit of the matrix model and 
its physical interpretation. To understand the effective geometry on Ai 2n , consider again a 
(test-) particle on A4 2n , modeled by some additional scalar field p (this could be e.g. su{k) 
components of (f) 1 ). The kinetic term due to the matrix model must have the form 

s[<p] = -^r[x^^[x^^7 ?a6 = -Tr([x^v 5 ][x^^^ + [0^^[^,^^) 

d 2n x p(x)G^(x)d„ipd u <p (31) 



(2 



7T 



where 

9itv(x) = rj^ + d^d^Sij = d^d^Tfab (32) 

GT(x) = e^\x)6 vu \x)g^ u ,{x) (33) 

= \0^\ 1/2 = \G^\ gixv (x)^. (34) 

Here g^x) is the metric induced onX 2n C M. D via pull-back of r] a b on R D . Now g^x) is 
no longer flat in general. So far, the kinetic term does not quite have the correct covariant 
form. This can be achieved by a suitable rescaling of G^ u (x): generalizing the corresponding 
quantities in @ and ([LTD, we define 

& v {x) = e- a G^(x) 
pG^ = \G^\ l ' 2 G^{x) (35) 

e- {n - l)a = |GV| 1/4 Mz)r* =pWi)l _i 

\G^\ = \e-lpt\g» v {x)\^ (36) 

Then the action (|3"T| ) has the correct covariant form 

S[<p] = [ d 2n x IG^G^Wd^d^. (37) 



7T 



Therefore the kinetic term on Aif 1 is governed by the metric G^ iv {x) 1 which has almost the 
same form as ( |TT1) except that the constant background metric g^ u is now replaced by the 
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induced metric g^, v {x) on A4 2n c M, D . The matrix model action ( p3j ) can be written in the 
semi-classical limit as 

S YM = -Tr[X a ,X b ][X a ' ,X b '] Vaa/Vw ~ ^ / (38) 

where 

477(2/) = G^(x)^(x) = (?» + 8^45^)6^' 6™'^, + dstfdstfSw) 

[X a ,X b ][X a ',X b ']r laa ,r ]w (39) 



generalizes ( jl6| ). 

There are 2 interesting special cases. For 4-dimensional NC spaces, we have 

\G^(x)\ = \g^(x)\, 2n = 4 (40) 

which means that the Poisson tensor 9^ u does not enter the Riemannian volume at all. 
This provides a very interesting mechanism for "stabilizing flat space" , and may hold the 
key for the cosmological constant problem as discussed below. In the case of 2-dimensional 
NC spaces, (|35|) has no solutionF], so that the action cannot be written in standard form at 
all. This will be discussed elsewhere. 

The emergence of such noncommutative vacua is very compelling in closely related 
(Euclidean) matrix models admitting compact NC branes as vacua [^,|30|], and supported 
by a considerable body of analytical and numerical work at least in 2 dimensions, including 
|H| |3"T|-|3"3]] and references therein. In higher dimensions, it may be necessary to consider 
supersymmetric matrix models as discussed below, cf. [34|. 



Relation with string theory. In string theory, a somewhat related situation occurs 
in the context of D-branes in a nontrivial B-field background. This leads to an effective 
description in terms of NC Yang- Mills theory on a noncommutative D-brane with Poisson 
structure Q^ v inherited from the B field, see e.g. [[27| and references therein. This effective 



gauge theory is governed by the open string metric which is strongly reminiscent of 
G^i^x) (apart from the density factor), while g^ v {x) corresponds to the closed string metric 
(more precisely its pull-back on the brane). Most of these results are restricted to the case 
of constant 9^ v and slowly varying fields, while the case of general NC curved branes has 
received only limited attention, notably [[35]. |36 . 



However, the results of the present paper should be compared more properly with 
previous work on string-theoretical matrix models such as the IKKT model |T2[ . NC branes 
have indeed been studied in considerable detail in this context, and it is well-known that the 



5 1 would like to thank A. Much for related discussions 
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matrix models can be interpreted as NC gauge theory on the brane. However, this has been 
worked out only for NC branes with a high degree of symmetry, such as fuzzy spaces (see 
e.g. or other special branes satisfying a BPS condition @ @, [23], [§] . The 



role of the effective metric (p3|) is well-known in these cases, and evidence for the existence 



of gravitons on the branes has been obtained [26]. For generic NC branes in matrix models, 



the effective metric G^ u (x) and its role in the effective field theory on branes has not been 
elaborated previously, to the best knowledge of the author. Moreover, it is essential to 
note that the would-be U(l) gauge field on the brane is absorbed in G^ v (x), leading to 
a dynamical emergent gravity. Therefore the present approach could be seen as a novel 
way of obtaining gravity from string-theoretical matrix models, avoiding the conventional 
picture of string compactification. 

In this context, it is worth recalling the relation between the semi-classical action (|3§|) 
and the Dirac-Born-Infeld action for 9~„ := B^ + F^ which governs the dynamics of branes 
in string theory |[27| . The action (|21|) arises from the DBI action at leading "nontrivial" 
order, 

yfdet( gtlu + 0rf) ~ p(x) (1 + 2 v {x) + ...) (41) 
omitting all constants, cf. p9|. 



Equation of motion for test particle ip. The covariant e.o.m. for <p obtained from 
the semi-classical action ([37]) is 

A d p = (G^d,d u - f><9> = . (42) 

On the other hand, starting from the matrix model fl3"T| ) we obtain the e.o.m. for the same 
scalar field f as 

= [X a , [X b , p]]Vab = [X», [X", (p])^ + [4>\ <p)} Sij 

= i[x», e^d v <p]n^ + i[4>\ e^d^dM 5 iS 

= -e^d p {e^d v <p)n^ - e^drfdrip^d^dvtp) <y - 

-O^in^ + Sg^dpd^ 

e.o = m. _Q mdp Q^ ( 4 3) 

The last equality holds for on-shell geometries defined by (f||), and Sg^ = d^ % d v <^ 5^. 
Comparing with the covariant form (§21), it follows that 



G^d^p = = A^, (44) 

for on-shell geometries, which implies the "harmonic gauge" 

f> c -°= m - 0. (45) 

This holds only in the preferred x M coordinates defined by the matrix model; a direct 
derivation based on (|48|) is given in Appendix A. 
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Equation of motion for X a . The same argument as above gives the equations of motion 
for the embedding functions <p l in the matrix model (|23|) , 



Aeft = (46) 

and similarly for x M ~ X^, 

Aqx* = 0. (47) 

This reflects the freedom of choosing the separation of X a = (X M , tp 1 ) into coordinates and 
scalar fields. In particular, on-shell geometries ( |48|) imply harmonic coordinates, which in 
General Relativity [40| would be interpreted as gauge condition. We will now derive an 



equivalent but more useful form of (|47|) in terms of the "tangential" 9 J; (x 



Equation of motion for 9 l{x). Reconsider the e.o.m. for the tangential components 
X M from the matrix model (|2Tj|): 



o = [x\ [x\ x b '}} nw = [x», \x\ x*\\nrt + [<t>\ [x v , n 5, : 



9^8 P 9^ Vwi > ~ 9^8^8,(9^8^ 8^ 



= -9" p d p 9^(r, m + Sg m ) - 9^9™d p 8g m 

= -9 w 'G^\x)8 p 9- v ^-9 v W ( 'd p g m (48) 

since d p 8g m (x) = 8 p g m (x), i.e. 

G pr} {x)d p 9~^ = 9^d p g^{x) = J u . (49) 

These are essentially Maxwell equations coupled to an external current J u , which depends 
on the matter field <p. As shown in Appendix B, this can be written in covariant form as 



G™{x) V 7 (e CT C) = e" CT 0™ 8^{x) (50) 



Here V denotes the Levi-Civita connection with respect to the effective metric G^ 
which is no longer unimodular in general. This has the same form as ( [15]) , and can be 
rewritten as 

V 7 0- x = -G m {x)9^8 1 a + e~ 2CT G, u 9^ d 7 r]{x) (51) 

The derivation in Appendix B assumes that the embedding functions (ff also satisfy their 
e.o.m. (|46|). It can also be derived directly from the semi-classical action (|38|): 

Semi-classical derivation of e.o.m. for 9~^(x). Starting from fl38|), we can derive the 
covariant e.o.m. of the matrix model using 

W£ = V^A, - V„<LV (52) 
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This gives 

SSym = 2 J d 2n x (s V (x) ^det Q-l + rj(x) 2 ^_L_ det 9^(9^89^)) 
= J d 2n xp(g, v 9^59™'g p>v , + g^'9 w '8g^, + rj(x) (9^89^) 
= J d 2n xp(G^9-lG^56^ + CTSg^ + v(x) (9^89^] 
2 / d 2n x yf\G](G r >i"G vp e' 7 9~^V p 8A ri - e^r/^V ' p 8A n 

+2 J d 2n x \[\G~\G» l 'd^dM Sij (53) 



using 

P 

which follows from (|35|). Noting that 

j2n 



l^l 172 e- CT (54) 



rf in x V|C| = (55) 



and VG = we obtain 



5S = -2 J d 2n x^J\d\5A v {G^G^V p ^9^) - V p (e^ V 0™)) + Stf 8^ (^ft] drf) 
= -2 J d 2n x^J\&\ [sA^G^G^V^e^-D - ^-^d^G^e-*^)) + <50%A 6S 
= -lJ d 2n x ^\ (^A^G^&PV p (e°9; t l) - e'^ d pV ) + % A^) 



using ( |180p and (Q) in the last steps. This gives precisely the equations of motion (^0|) 
and 



Formal considerations. From a more formal point of view, we have the following struc- 
tures: The submanifold Ai 2n C M D carries an embedding metric g, and a preferred frame 
e fi = Q& v Q v r*j [X^, .] which encodes the noncommutative structure. The effective metric 
G pv (H) on T*M is defined by 

(P, P') G := P, P' E T*M , (56) 

where (d p , d v ) g = g pu (x) and * : T* M. — > TM. is the canonical map defined by the Poisson 
structure 9 pv . 

Notice the unusual role of the indices. This makes sense here, because the frame e p is 
given in terms of the antisymmetric Poisson structure 9 pv in the preferred coordinates x p . 
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There is no distinction between "Lorentz" and "covariant" indices here, and neither local 
Lorentz nor general coordinate transformations are allowed a priori. One could proceed 
to introduce differential forms in terms of one-forms 9 a , a = 1,...,D through [^X 6 ] = 
0, 9 a 9 b = —9 b 9 a . The exterior differential of functions is then defined in terms of a "special" 
one- form 9, 

df=[9J], 9 = X a Vab 9 b . (57) 

This is similar^ to the formalism in |4l| , |42"|| ; however the calculus is D-dimensional, similar 
to the case of the fuzzy sphere El. The scalar action can then be written as 



S[tp] ~ / d 2n xp(d(f),d(f)} (58) 
where (9 a , 9 b ) = i] ab . Similarly, the semi-classical form of the matrix model is 



Sym ~ 777 
(Z7T 



^— / d 2n xp(9 A9,9 A9). (59) 

K) n J 



We can also write 



df = [9,f] =e"(/)0 M , 

K = V,u9 u + d^9 2n+ i (60) 

where 9^ is in some sense dual to e M . This should illuminate the relation and difference 
to [[|T[]. These considerations will be pursued further elsewhere. 



2.2 Nonabelian gauge fields 

Now consider backgrounds of the form 

Y a = X a ® ln + A a a ® X a (61) 

where A" are generators of su(n). According to [ffl, the U(l) sector (i.e. the components 
proportional to l n ) is absorbed in the geometrical degrees of freedom defined by X a , and 
the discussion of the previous sections applies without change. On the other hand, the 
su(n) components A% behave as nonabelian gauge fields, and similarly the transversal 
su(n) components (p l a in 

<^ = <^ <8> 1„ + ® A a (62) 

are nonabelian scalars from the brane point of view. The <j) l a then propagate in the back- 
ground geometry as discussed above. If some of the cf> l a develop a nontrivial vev, they 
might be viewed as part of the geometry. 

6 However the frame and metric here have a specific form in terms of Q^ LV , unlike in |4l| . 
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It was shown in 0] that the effective action for nonabelian gauge fields A£ due to the 
4-dimensional matrix model (H) in the semi-classical limit is 

S YM [A] ~ J d 4 x p{x)tr (G^'G W 'F^ F„ v ) + 2 J rj{x) trF A F 

= J d 4 x |G^| 1/2 e CT G^'G uv 'trF^ F^ + 2 J r](x) trF A F . (63) 

This is the Yang-Mills action for a nonabelian gauge fields coupled to the effective metric 
G^, apart from the "would-be topological term" and the density factor e° ' . The latter 
could be interpreted as varying bare gauge coupling "constant" 

9 2 ym = (64) 

introducing an overall coupling constant to the matrix model ([I]). In order to be phys- 
ically acceptable, it is probably required that a is slowly varying. Indeed, a kinetic term 



for the "dilaton" p resp. a is induced in the quantum effective action except in the 



case of unbroken N = A supersymmetry. This might ensure that a is nearly constant. 

Due to the strong constraints of gauge invariance, we expect that ( |63"D applies without 
change to the case of non-trivially embedded 4-dimensional branes in M D ; however this re- 
mains to be shown. Note that rj ~ e u due to (jS8|) , hence the two terms in (|63D have roughly 
the same coefficients. This changes for higher- dimensional branes, where the "would-be 
topological term" J rj(x) trF A F will be replaced by a different term which could be deter- 
mined along the lines in ||]. Before relating this e.g. to the strong CP problem one would 
first have to identify more realistic models, elaborate the symmetry breaking etc.. 



2.3 Fermions 

Then the most obvious (perhaps the only reasonable) action for a spinor which can be 
written down in the matrix model frameworkF] is 

S = (27r) 2 "Tr^ 7a [X a ,vl/] „ J d 2n x p(x) *»( 7a4 + l2n+i d^ i )B^{x)d v ^ 

= J d 2n xp(x)m%6^(x)d u y (65) 

where 7a defines the D-dimensional Euclidean Clifford algebra, and 

% = Ifi + l2n+id^(p l (66) 

satisfies the Clifford algebra associated with the embedding metric g^x) on A4, 

{%, %} = 2^ u + 2d^ dvtfSij = 2g^{x) . (67) 

In particular, fermions should also be in the adjoint, otherwise they cannot acquire a kinetic term. 
This does not rule out its applicability in particle physics, see e.g. Q. 



13 



This is indeed the appropriate coupling of a spinor to the background geometry with metric 
G^v (up to rescaling), albeit with a non-standard spin connection which vanishes in the 
x M coordinates. This nicely generalizes at the classical level the analysis in NTT], where 



this action was shown to provide a reasonable coupling of fermions to emergent gravity 



for flat g^ u . At the quantum level, it was shown in |l Lj that the Einstein-Hilbert term is 
indeed induced (along with a Dilaton-like term for a), for flat g^ u and on-shell geometries. 
It remains to be verified whether this generalizes to the case of non-trivially embedded 
branes. This is expected to be the case since it does give the correct Dirac operator e.g. in 



the case of the fuzzy sphere EH] or for Sff x S* 2 



A? 



29 



Given the above Dirac operator, one could also consider the associated spectral action 
in the sense of ||46| . It is an open problem in that context how to quantize gravity, more 
precisely how to integrate over the various geometries. The present framework suggests a 
simple answer: The Dirac operator should have the form as given in (|65|) , and the integral 
over the geometries should be realized as integral over the matrices X a with measure defined 
by the bosonic matrix model, d/i(X a ) = e~ SyM ^ xa ^ (p3|). Nevertheless, this is not entirely 
equivalent to the present matrix model framework: The spectral action is based on the 
dependence of the spectrum as a function of the cutoff, while in the N = 4 considered 
here such a cutoff should not be required. 



2.4 Compactification of branes 

Consider a 2n-dimensional NC brane Aig n C M 10 . In order to obtain a 4-dimensional 
space at low energies, we assume that this higher- dimensional brane has compact extra 
dimensions, for example 

Mf ~ M.q x Kg . (68) 

If K is "small" enough, this looks like M.\ at low energies, as in standard compactification 
scenarios. Particularly natural examples would be ~ Mq x or JHq ~ Mg x Sj^ x 
Sjf, where S% denotes the fuzzy sphere. Such extra-dimensional fuzzy spaces can indeed 
be embedded naturally in the matrix models considered here (possibly upon adding soft 



SUSY-breaking terms) [|j|,|3T| , or alternatively they can arise spontaneously from the scalar 
fields from the 4-dimensional point of view [f44| , |4"7|l . These 2 points of view are essentially 
equivalent. 

Let us count degrees of freedom for the effective metric. For a 2n - dimensional NC 
brane, 9^ v resp. 6~ l = dA has 2n — 2 physical (on-shell) plus one off-shell degrees of 
freedom, after gauge fixing. Upon compactification on Kg, the components A\ tangential to 
Kg become massive, leaving only 2 massless d.o.f. from a 4-dimensional point of view. The 
embedding of M 2n C M 10 defined by (ft 1 provides 10 — 2n additional degrees of freedom. They 
are absorbed in the effective metric and governed by the quantum effective action. From the 
4-dimensional point of view, this will lead to an effective "brane tension" on Aig depending 
on the moduli of the compactification (e.g. the radius) as indicated below. Those are likely 
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to become "off-shell" d.o.f. which enlarge the class of effective 4- dimensional metrics, as 
desired. Therefore one can expect to recover most of the 2 on-shell plus 4 off-shell d.o.f. of 
the 4-dimensional metric in General Relativity. All this requires a more detailed analysis. 



Example 1: The fuzzy sphere Sfj. The fuzzy sphere Sfj [43| is a natural realization 
of this framework, being realized in terms of an embedding S 2 cK 3 . Consider our matrix 
model in D = 3, with the configuration 

X a = JL X {N) ' a , a = 1,2, 3 
c N 

[X a ,X b ] = i6 N e abc —, X a X a '5 aal = ^-(N 2 -l) = r 2 . (69) 

r c N 4 

Here r is an arbitrary radius, and \( N >> a denotes the generators of the N- dimensional 
irreducible representation of SU(2), 

4 = ^ 2 -i) (70) 

0n = —■ (71) 

C N 

Even though ( |69|) is not a solution of the basic matrix model ([!]), it makes nevertheless 
sense to consider such configurations since the induced gravity action is not yet taken into 
account. Moreover, it becomes a solution once a mass term and/or a cubic term is added 
to the action, e.g. 

SVm + 5 corr = (2n)Tr([X a ,X b }[X a ',X b '}5 aa ,5 w + m 2 X a X a ' 5 aa , + 7 X a X b X c e ab ^j . (72) 

Such terms might be induced in the quantum effective action, possibly after SSB. The 
general discussion of Section ^[1] applies as follows: consider e.g. some neighborhood of the 
north pole x 3 « r, i 1 « x 2 <C r of S 2 . Then we separate the coordinates as in (|24"D in 2 
tangential ones and one scalar "embedding" function, 

X a = (A^,0), n=l,2 (73) 



where <fi = 4>{X^) m for small X l ,X 2 . Indeed it is not difficult to write = X 3 in (|69|) 
as a function of X 1 ^ 2 , 
The Poisson tensor e.g. c 
the symplectic structure 



as a function of X 1 ^ 2 , in a suitable domain. Hence S% is a NC brane embedded in R 3 . 
The Poisson tensor e.g. at the north pole is 6 12 = 6n, and S% is a quantization of S 2 with 



u S 2 = 9 N l - e abc x a dx b dx c , (74) 

where x a is the semi-classical limit of X a . It satisfies the semi-classical quantization con- 
dition 

>s 2 Js 2 



/ u S 2 = / d 2 x p = 9 N X Anr 2 = Atcc n (75) 
Js 2 Js 2 
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consistent with (ffD[), where p = O^ 1 . Therefore can be considered as a compactification 
of the D = 2 Moyal-Weyl plane. The embedding metric g^uix) is the round metric for a 
sphere S 2 with radius r, and G^ v = 2 N g^ v . 



Example 2: Rg x S N . Now consider a configuration Aig = M.gX Sfj. This can be realized 
in the D-dimensional matrix model for D > 7 

X» = X^, n = 0,1,2,3 

ft = i = 1,2,3 (76) 

cn 

where X^ are the generators of JR| ( |105| ). This should be interpreted as a 6-dimensional 
NC space, which for small r looks like M^. Such configurations can lead to interesting low- 



energy gauge groups and zero modes in the nonabelian discussed in [Tj.nj . Similar 

configurations were discussed previously in the IKKT model see also J48|. The radius 
of the fuzzy spheres will be dynamical r = r(X M ) ~ r(x M ), determined by the effective 
action. Inserting this configuration in the action and recalling ([3^), we obtain 

S YM = (27r) 3 Tr[X a ,X b }[X a ',X b '}6 aa/ d bb ,^4[ d 4 xu s2 p^\x) v ®{x) 

JR 4 x5 2 

= 8ttN f d 4 xp {4 \x)7]^(x) (77) 
using (|75|) and 

^\x) ~ [X a ,X b }[X a ',X b '}S aa/ S w ~ ^\x) + 2G""^r(x)d i ,r(x) + 2^, (78) 

C N 

cf. (H). Here 

V W( x) = G ^g^ x) (79 ) 

involves only the 4-dimensional metric. This leads to an effective potential V(r) for the 
radius r(x), which will receive additional contributions from further terms such as ([72] ) and 
from the induced gravitational action. For example, consider the 6-dimensional "would-be 
cosmological constant" term in fl96), which using ( |36"D can be written as 




flU 



Anctf \ d 4 xr{x)p{xfl 2 \g ilv \ x l 4 A % (80) 



where p(x) and g^x) in the last line are 4-dimensional quantities. Somewhat surprisingly, 
this term now depends on p(x), unlike in the case of 4D branes. However, this expression 
should be taken with much caution, because the IR condition ( |90|) for the applicability of 
the semi-classical expressions ( p6|) may not be appropriate for the compact dimensions. In 
that case, it might be more appropriate to use a 4-dimensional description rather than the 
6-dimensional metric. In any case, this will contribute to the effective potential for V(r), 
but a more detailed analysis is required. 
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Example 3: M.g x Sff L x Sfj R . A generalization of the above configuration which can be 
realized in the 10- dimensional IKKT matrix model is 



X» = X», fi = 0,1, 2, 3 

ft = Hx^'\ i = 1,2,3, 
cl 

ft = ^ X (N R ),i i = 4)5j6 
CR 



which should be interpreted as a 8-dimensional NC space. The effective 4-dimensional 
action now involves 2 parameters tl,tr, which will be governed by an effective potential 
V{tl, tr). This should provide sufficient structure to obtain interesting solutions from the 



particle physics point of view; see also e.g. [29,31] and references therein for related work. 



2.5 Departures from General Relativity: preferred scales and 
coordinates 

There are several features of the model under consideration which differ radically from 
the conventional picture of General Relativity. We focus on the case of 4-dimensional NC 
branes for simplicity. 

First, recall that there are preferred coordinates in the model, given by the covariant 
coordinates x M . In those coordinates, the background metric is explicitly constant, g^ v = S^ u 
resp. = rj^ in the D = 4 case, and the preferred frame is given by the antisymmetric 
Poisson tensor, e M = 9^ v d v . This is physically not very significant a priori, but it simplifies 
the issue of gauge fixing. We recall that in the preferred x M coordinates, the on-shell 
condition for amounts to ( |45|) 

f i* = , (82) 
which would be interpreted as gauge choice in General Relativity. 

A more significant feature of the matrix model is the presence of the scalar density 
p = (det Ouv) 1 ^ 2 (P^D, which defines the scale of noncommutativity 

p = A% c = L~ N 4 C (83) 

and provides the symplectic measure (27r) 2 Tr/ ~ f dxpf. Such a structure does not 
exist in the commutative framework. This leads to an analog of the Bohr-Sommerfeld 
quantization condition, 

Vol e = {2ix) 2 N (84) 

where Volg denotes the volume measured in units of L^c, and Af the dimension of the 
corresponding Hilbert (sub)space. This means that the volume is quantized in integer 
multiples of L A NC , so that NC branes are automatically "large" for large Af. This is already 
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a hint that NC spaces like to be flat, which is very interesting in connection with the 
cosmological constant problem. 



There is another scale in the model determined by the embedding metric resp. the 

Lj = A; 4 = | | 1/2 (85) 



effective metric G^, 



which we could set to 1 thereby fixing the units; recall that \G^ IU \ = ( ffDD for general 
4-dimensional branes in ~R D . The ratio of these scales defines the dimensionless scalar 
function 

\e,jv 2 At 



\G^ A 4 g K } 



using (|3q) . We can relate this with the Riemannian volume of (-M$, G^ v ) measured by G^ u , 

Vol d = {2n) 2 ATe a . (87) 

The "dilaton" e CT will be determined dynamically by the model resp. the background 
under consideration. For example, in matrix models for fuzzy spheres it depends on the 
coefficient of additional (soft SUSY-breaking) terms such as Tr e a b c X a X b X c , see also the 
related discussion in ||7). Note that e~ a also gives the scale of 77, 



\G(j.u 


1/2 


\0fiu 


1/2 



V 

at least for simple 4-dimensional configurations. This may be a significant large dimension- 
less number. 

In the context of quantization, we will encounter 2 additional scales A4 ^> Ai in the 
10-dimensional version of the model, where A4 is the scale of A = 4 SUSY breaking which 
is argued to coincide with the Planck scale Api below, and A x is the scale of N = 1 SUSY 
breaking. These should also be dynamical scales. We will furthermore argue that A^c > A 4 
simplifies the semi-classical analysis, however this is not essential; it seems actually plausible 
that N = 4 SUSY is broken by the NC background, so that A^c = A 4 . In summary, we 
expect 3 a priori distinct physical scales 

A NC > A 4 = /p/ > Ai (89) 

in addition to the dimensionless number e a in the model. 



2.6 Quantization and induced gravity 

Now consider the quantization of our matrix model, which can contain scalar fields (such 
as e.g. arising from extra dimensions), fermions, the "would-be U(l) gauge field" which 
is absorbed in 6"^ (2), and possibly nonabelian gauge fields. In principle, the quantization 
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is defined in terms of an integral over all matrices. This is expected to be well-defined 



at least in the case of the IKKT model, which leads to N = 4 SUSY on M.g [13 



Some modifications such as soft SUSY breaking terms may also be allowed. Note that 
this quantization implies an integration over all geometries of the NC branes embedded in 
M 10 , via the quantization of the embedding functions (jf as well as X M resp. 6 >J,v {x). In 
particular, (emergent) gravity is also quantized. 

To obtain a qualitative understanding of the model at the quantum level, we can take 
advantage of the above semi-classical form of the action in terms of conventional field theory 
coupled to G^ v . Then the low-energy effective action at one loop can be extracted from 
standard results of ordinary quantum field theory on curves spaces. As shown in |fL0[, jn|, 
this is indeed justified (based on a comparison with a fully NC computation and UV/IR 
mixing) provided there exists an effective UV-cutoff A, and the following IR regime |10[ is 
respected 

pA < A 2 NC and A < A NC ■ (90) 

These conditions ensure that the effects of noncommutativity are mild even in the loops, so 
that the phase factors in non-planar diagrams are small and are well approximated be the 
Poisson structure. This reflects the fact that emergent NC gravity is an IR phenomenon. 
A violation of e.g. A < A^c is acceptable, but implies corrections^ to the effective action 
(f)6|) given below, some of which have been discussed in flTol . Such a cutoff is realized in 
the N = 4 supersymmetric version of the model, assuming that N = 4 SUSY is broken 
at A = A 4 from now on. This is essential, because no bare term in the action is available 
which could cancel the induced (gravitational) action discussed below. We will furthermore 
assume that some smaller supersymmetry survives down to a much lower energy scale Ai, 
below which no supersymmetry survives. These are reasonable assumptions, which appear 
to be necessary for the proposed framework to be physically viable. Note that these scales 
are measured using the physical metric G^ v . 

The results of the one-loop computation of fields coupled to the background metric 
can be obtained conveniently using the Seeley- de Witt coefficients of the corresponding 
heat kernel. The essential features are illustrated by the quantization of scalar fields. Hence 
consider the effective action obtained by integrating out the scalars, which in the Euclidean 
case is 

e -W= / rf$e" 5 ^. (91) 



Since we are mainly interested in the induced gravitational action here, it is sufficient to 
consider the case of non- interacting scalar fields coupled to the metric G^ Vl where 

r^^Trlog^Ae (92) 



8 if this condition is violated, a more refined analysis of NC corrections is required, cf. . It turns out 
that the apparently-quartic divergent term J Af \f~£} actually becomes milder, being a difference between 
quadratically-divergent planar and non-planar diagrams. A similar comment applies to the J A|i?[G] term. 
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assuming Euclidean signature for simplicity. Here A^, is the Laplacian of a scalar field on 
the Riemannian manifold (M, G ab (y)) with action fllQl) . The UV cutoff A is incorporated 
using the Schwinger parametrization 



Trflog^A^ _ log^A 



o " 
00 da 



= _Tr / — [ e -°5 A G _ e -4M e -^2 



a: 



(93) 



Now we can use the heat kernel expansion, 



m>0 



(94) 



The a m (x,A ( g) are known as Seeley-de Witt (or Duhamel) coefficients, which for scalar 
fields with action fllOf) are given by [49 



f 



a (x) 
a 2 (x) 



1 /l 



(47r) n V6 



R[G]), 



(47r) n 360 



(95) 



The full effective action is complicated by the fact that there are several types of fields 
in the model including scalars, fermions, and gauge fields. Because they all couple to the 
same effective metric up to possibly density factors^, they will all induce essentially the 
same type of gravitational action, with an additional kinetic term for the "dilaton" a due 



to the fermions |TT| and gauge fields resp. gravitons. Therefore one obtains the following 
type of induced gravitational action: 



Fl-loop[G] 



(47r) n 



d 2n x 



G 



fJtU I 



coAf -c 2 R[G) Af- 2 + 



(96) 



where c m are model-dependent constants, omitting dilaton-like terms. This allows already 
to draw some qualitative conclusions, focusing on 2n = 4 -dimensional NC branes. More 
detailed computations should be performed elsewhere. 



Note that we associated different scales Ai resp. A 4 to the different terms in (|96|), which 
arise as follows. It is well-known that the coefficient of the leading "would-be cosmological 

constant" term J d 4 x \J \ G^ v \ is determined by the scale Ai for N = 1 SUSY breaking. In 
contrast, the coefficient of the induced Einstein-Hilbert term in emergent NC gravity is 



9 This was shown in for gauge fields and in pi] for fermions. Similar results are expected for the 
gravitons (i.e. the would-be U(l) gauge field) due to supersymmetry, at least in the case of N — 4 SUSY. 
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determined by the scale A4 where N = A SUSY is broken. This reflects the well-known 
fact that UV/IR mixing in NC gauge theory persists even in SUSY gauge theory |50"|, jSlfl , 
except in the N = 4 case. Since ri_j oop [G] is nothing but a re-interpretation of the UV/IR 
mixing terms in NC gauge theory 0,[lO|, it follows that the induced term R[G] has a cutoff 
given by the scale of iV = 4 SUSY breaking; this is discussed in [11] from the point of view 
of gravity. Because there is no bare gravity action, it follows that the effective Newton 
constant resp. the Planck scale in emergent gravity is given by 

& = i ~ Al . (97) 

This also suggests what happens in models without a finite effective cutoff A 4 : Then G — > 0, 
hence the the induced gravitational action becomes a constraint and there is no more back- 
reaction of matter to the geometry. However there might still be interesting scaling limits. 



2.7 Effective action and the (ir)relevance of the cosmological con- 
stant term 

Consider again the case of 4-dimensional NC branes embedded in IR D . The full semi- 
classical effective action of the matrix model at one loop is given by 

$eff — Sym + Si-i 00 p (98) 
where Si-i oop = IV^fC] + and 

Sym = J d 4 x^\G~^\ (e^G^'G^trF^F^ + G^'trd^d^Sij) 

, ~vT~5 f d * x ( A PV + M x ) trFAF) (99) 
(2tt)V J 



+ 



including nonabelian gauge fields and the nonabelian components 4> l = <p l a X a of the scalar 
fields. Fermionic terms are omitted. We introduced an explicit coupling constant g, which 
does not enter the induced gravitational action since it can be absorbed in the fields. The 
bare YM coupling constant is given by g Y M = Q 1 ^ (0)) which receives the standard 
quantum corrections, and might play a role similar to a GUT coupling. The one-loop 
induced gravitational term Ti_i oop [G] for 4-dimensional NC branes is fl96|) 



ri_ JO op[G] ~ J d 4 x (At + R[G]Al) (100) 

where = {g^ul using (|40|). The first term is therefore simply the invariant volume of 

the embedding metric. 

Consider the geometric equations of motion. It is well-known that 

S J d 4 x R[G] = J d 4 x yJ\G~\ (^R[G]G» V - R^fi^SG^ (101) 
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while the variation of the "would-be cosmological term" is 



5 J d 4 x ^J\G~\ =\J d ' x \[\G~\G^ SG^ = \j d A x JggT 5g^ v , (102) 

using ( f40"D in the case of 4-dimensional NC branes. This vanishes identically in the case 
D = 4 due to (|TT|) . The variation of the bare gravitational action 5 J d 4 x4pr] was worked 
out in (1561). 



In General Relativity, ( |101| ) and (|102| ) imply the Einstein equations for vacuum. The 
essential difference here is that the metric G^ v is constrained, and the fluctuations do not 
span the space of symmetric 4x4 matrices. Therefore we do not simply obtain the Einstein 
equations; this is seen most strikingly for the cosmological constant term as discussed below. 
However, note that Ricci-flat spaces which can be realized in terms of our G^ v certainly 

satisfy 5 J d 4 x J \G\ R[G] = 0. This supports the physical viability of this framework. The 
correct e.o.m. which follow from the effective action are complicated here by the presence 
of possible dilaton-like terms at one loop, and will be derived elsewhere. They will in 
particular modify fl5Pp, which takes into account the bare action only. 

Let us briefly discuss the geometrical degrees of freedom. We can decompose the vari- 
ations 5X a of the basic matrices into tangential and normal fluctuations w.r.t. the back- 
ground brane Ai. Using an orthogonal transformation if necessary, we can assume that 
5X^ G TA4 are tangential, and 5(f) 1 G TAi 1 - are normal to the brane at some given point. 
Consider first the tangential variations 5X^ = A^ 1 . They lead to variations of the Pois- 
son tensor 6^ u on the brane (which can be interpreted as diffeomorphism) , but they do 
not change the embedding metric, which is fixed in the matrix model (recall that e.g. 
g^ u = rj^y in the simplest case D — 4): 5j L g iW = 0. Therefore these tangential fluctuations 
imply nontrivial|3 physical fluctuations of the effective metric 5a G^ v ~ (G9F + FOG)^ 
( p.!5| ) corresponding to the 2 on-shell graviton helicities plus one off-shell deformation; cf. 

the discussion in Section |3|. In particular, the term J d 4 x \J\G\ = J d A x \/\g~\ in ([51]) is 
independent of these tangential degrees of freedom. This provides (part of) a mechanism 
for avoiding the cosmological constant problem. 

Now consider the normal fluctuations 5(J) 1 of the brane embedding, which in general 
imply nontrivial physical fluctuations of the effective metric. The corresponding variation 
of the "would-be cosmological term" is 



5 J d A x = J d A x ^W\g^5g, u = 2 J d A x 'd^djfi % 

= -2 J dtxy/lgiSftAgfiSij (103) 
using partial integration, where A g is the covariant Laplacian corresponding to the metric 



10 



except for gauge transformations resp. symplectomorphisms = [/, X^]. 
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g^y. This vanishes if the (j) 1 satisfy the constraint 



A ff 0* = O, (104) 

which is similar to ( [46]) except that the metric is now g^ u . Flat embeddings do satisfy this 
condition. Therefore flat space is a solution at the quantum level, even in the presence of 
this "would-be cosmological constant" term; the same applies to any surfaces embedded 
in M D which satisfies (|104|) . This can also be seen from the fact that there is no tadpole 
contribution at one loop in NC gauge theory [[Hj . This is in stark contrast to General Rela- 
tivity, where the term J <i 4 x\/(j A 4 corresponds to a huge cosmological constant, requiring 
unreasonable fine-tuning. Together with the observations of the previous paragraph, we 
obtain strong evidence here that the cosmological constant problem is resolved or at least 
much milder in the present context. This is a robust mechanism, rooted in the fact that 
the metric is not a fundamental degree of freedom but emerges from the matrix model. 

To obtain a more complete understanding of the cosmological constant issue in emergent 
NC gravity, a more complete analysis is required, as for related claims in the literature [[)2[ 



At present, the only known solution for the full effective action (j98|) is flat Moyal-Weyl 
space. Nevertheless, the fact that this is a solution without fine-tuning a cosmological 
constant is very remarkable. Moreover, since the Einstein-Hilbert term contains two explicit 
derivatives, the bare action together with the "would-be cosmological constant" term will 
govern the extreme IR (cosmological) scale which should indeed be flat, while the induced 
E-H action will determine the gravitational fields due to localized (point) masses. This 
would be a very satisfactory picture. 



3 Linearized metric and gravitational waves 

Moyal-Weyl case. A particular solution of the e.o.m. (|T3|) is given by the 4D Moyal- 
Weyl quantum plane. Its generators X^ satisfy 

[X* i ,X v \ = i6' a 'l, (105) 

where is a constant antisymmetric tensor. The effective geometry (]?]) for the Moyal- 
Weyl plane is indeed flat, given by 

gT = pgT 

p = \g^\ 1/4 = \e-J\ 1/2 ^ A% c . (106) 

In this section, lower-case g^ u resp. g^ v will denote the flat effective Moyal-Weyl metric 
rather than the embedding metric, and we will rise and lower indices using g^ v . First, we 
can choose coordinates where g^ v = (—1, 1, 1, 1), so that x° = ct corresponds to the time. 
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One can use the remaining 5*0(3, 1) (resp. £0(4) in the Euclidean case) to bring 9^ u into 
canonical form 

/ -1 \ 
a 

-a 
\ 1 J 



9^ = 6 



(107) 



The NC scale is then 



A 4 

V/VC' 



p = e- 2 a x 

and the original flat background metric is 

V»v = 8~lKv>r' v ' = P' 1 ^ 2 diag(l, a~ 2 , a~\ -1) = a diag(l, a~\ ec\ -1) . 
The bare action for this flat Moyal-Weyl background is 



S YM = Tre^e*" Vw'Vw = Tr Vflu g 



2 / d x a (a 



r 



(108) 
(109) 

(110) 



which vanishes in the case a = 1 where admits an enhanced SO(2, 1) x U(l) symmetry. 
In the Euclidean case, the action is positive definite. From now on, we assume that 



a 



for simplicity. It is also worth pointing out that we are not in the case of "space-like" 
noncommutativity, since 6^ v is non-degenerate. However, the problems of unitarity etc. 
discussed e.g. in |k| are expected to be benign in the present context due to the assumed 
N = 4 supersymmetry at the Planck scale. 



Deformations of the flat Moyal-Weyl plane. Consider small deformations of the flat 
Moyal-Weyl plane, 

i^F-r^x), (in) 

where A v are hermitian and can be interpreted as 27" (1) gauge fields on M 4 , with field strength 
= d fl A u — d fl A u + i[A^, A u }. The linearized metric G^ u is 

= g^-h^ (112) 

where 

ir = -g^' p^ev - 9^' sr - o^'e^sg^ + o(a 2 ) . (113) 

where 5g^ u = d^ l d v <^ Sij. Correspondingly, the inverse metric is 

G lJLV = g lJLV + h^ v + ... , (114) 

with 

V = Muwh^ = -g vv S v ' p F pli - g^r' p F pu - g^g^'^^'Sg^. (115) 
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This gives 

h = K v g» v = 20*" - rTSg^ . (116) 

We focus on the case of flat embeddings Sg^ = 0. Then W v = -gw' F p , p 9 pv -9^ F p , p fp + 
0(A 2 ) gives the linearized fluctuation resp. graviton in terms of the U(l) degrees of freedom. 
The linearized Ricci tensor for the unimodular metric G^ u resp. the traceless graviton 

V = V - fv^ is g iven b y 

Rfiu[G] = i f V d^i^, + ^ c^i^ + \g^d p d p F n J^ (117) 
in agreement with results of H, and 



R[G] = ^d p 9 r ' a F ria . (118) 



Now consider the equations of motion for the bare action (14), which in the present context 
amount to d p F pv = = d p d p F fJiU up to possibly corrections of order 9, i.e. the vacuum 
Maxwell equations for the flat metric ~g . As pointed out in ||, this implies that the 
vacuum geometries are Ricci-flat to leading nontrivial order, 

R ab [G) = + O(9 2 ), (119) 

while the general curvature tensor Rfj, vpv is first order in 9 and does not vanishP). This 
shows that the effective metric does contain the 2 physical degrees of freedom (helicities) 
of gravitational waves. It is quite remarkable that (fTT^) is obtained from the bare action, 
without invoking the mechanism of induced gravity in section |2.6| . Note that the cosmo- 
logical constant vanishes to this order. A generalization to non-trivially embedded branes 
remains to be elaborated. 



3.1 Newtonian Limit and relativistic corrections 

The Newtonian limit of General Relativity corresponds to static metric perturbations of 
the form 

ds * = _ c ^ 2 (l + ^) + dx 2 (l + O(I)) (120) 

where A^U(x) = 4iiGm(x) and m(x) denotes the mass density. Including the leading 
relativistic corrections, this takes the form 

ds * = _ c 2 dt 2^ + ^) + dx 2 (l - ^) . (121) 

11 while this is true generically, there may be particular momenta determined by Q^ v for which the 
corresponding "graviton" is pure gauge and hence R pV pa vanishes. This should be studied in more detail 
elsewhere. 
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It follows from the results of the previous section that ( |121|) is reproduced correctly in the 
case without matter, where the vacuum equations of motion amount to d c F c h = resp. 
R a b = 0. In the presence of matter, it was essentially shown in [[| that one can indeed 
obtain metrics of the form ( |120| ) for arbitrary static m(x). We will re-analyze this issue in 
more detail here. It will turn out that even though one can find a metric corresponding 
to ( |120|) in the case D = 4 without nontrivial embeddings, the relativistic corrections of 



General Relativity in ( |121| ) are not correctly reproduced in the presence of matter. In 
particular, it appears that the Schwarzschild solution is not correctly reproduced in this 
minimal framework. This is not a real problem for emergent NC gravity, since we concluded 
on different grounds above that the model with D = 10 and N = 4 SUSY is required for a 
consistent model at the quantum level. Therefore realistic solutions for point masses should 
be realized by nontrivially embedded branes. 



In order to reproduce (|120|) with G^ u , we have to find U{1) gauge fields on the 
Moyal-Weyl quantum plane with the desired h^ iu , which is ( |115| ) 



°fJ,U 



r\v ■ 



:i22) 



Choose coordinates where g^ u = (—1, 1, 1, 1) as discussed above, so that x° = ct corresponds 



to the time, and 9^ u has the form (|107|). Then 



I 





E x 


E 2 


E 3 


\ 




-E x 





B 3 


-B 2 






—E 2 


-B 3 





Bx 




V 


—E 3 


B 2 


-Bx 





/ 



(123) 



gives 



h 



p6 



( — 2E3 B 2 — E 2 —Bx + Ex \ 

B 2 — E 2 — 2S3 Bx + Ex 

—Bx + Ex — 2.B3 B 2 + E 2 

\ Bx + Ex B 2 + E 2 2E 3 ) 



(124) 



which is the most general metric fluctuation available. Let us denote its trace with 

h(x) = r u K„(x) = A9(E 3 - B 3 ) . 
The physical graviton is the traceless version, 



(125) 



:9^uh = p9 



( -(B 3 + E 3 ) B 2 -E 2 -Bx + Ex \ 

B 2 — E 2 —(-83 + E 3 ) Bx + Ex 

-Bx + Ex ~{B 3 + E 3 ) B 2 + E 2 

\ Bx + Ex B 2 + E 2 B 3 + E 3 ) 



(126) 



As shown above, R ilv \G\ = holds if E and B satisfy the Maxwell equations without 
source. We would like this to be the case where the mass density m(x) vanishes. 
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3.1.1 Static point charges 



To gain some intuition, we first consider a static point charges with electric and magnetic 
charge at the origin, and determine the corresponding effective geometry. Thus consider 
electromagnetic fields given by 



E; 



1 



Bi 



X i 



(127) 



Then the metric fluctuation (|126 ) is 



( ~{Qe + Qm)x3 (qn ~ qE)x 2 

-{mm-Qe)^\ o 

\ (q M + q E )xi 



-{q.M ~ 1e)x\ 




o \ 

(qM + qE)x\ 



-(qe + qM)x 3 (qM + qE)x 2 

(qM + qE)x 2 (qE + ?m)^3 / 



(128) 



which in the "extremal" case qM = qE ='■ q is 



h 



2qp6 



( 



\ 



-X 3 









-x 3 


X\ 



\ 

x\ 

-x 3 x 2 

x 2 x 3 ) 



(129) 



This can be brought into diagonal form using the diffeomorphism ^ = 2qp9(0, 0, 0, ^ 
which gives the linearized metric 



ti. 





f 


-x 3 








o \ 


2qp6 







-x 3 


















-x 3 







\ 











-x 3 J 



(130) 



This has indeed the form of ( 121 ) of a Ricci-flat metric for f ^ 0, with Newtonian potential 



U — x 3 ~ d 3 - which is harmonic away from the origin. However, the corresponding 



'mass" distribution 



mix) 



AU ~ d 3 5 {3) ( 



x) 



(131) 



is not positive. This is not what we want; it corresponds to an unphysical gravitational 
dipole rather than a point mass. Note however that there is no charged field in the model 
for this U(l), hence this is only a toy configuration which is not expected to play any 
physical role. Moreover, it is not expected to be a solution of the e.o.m. at the quantum 
level. 

This result is easy to understand: Since the electromagnetic field of a localized charge 
distribution decays as \, the corresponding gravitational field also decays like \ at the 
linearized level. The correct U(r) ~ i gravitational potential for a point mass can be 
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recovered either at the cost of violating relativistic corrections as elaborated next, or - 
presumably - by a nontrivial deformation of the brane embedding due to the point mass, 
governed by the induced gravitational action. The trace-£/(l) modes under consideration 
here are to be interpreted as gravitational waves. 



3.1.2 General mass distributions 



Now consider more generally the static case 

h 0i = 0, <9 /v = * 

which amounts to 

B 2 = E 2 , Bi = E x . 
The metric fluctuation (|126|) is then 



(132) 
(133) 



( -(B 3 + E 3 ) 



V o 



\ 

[B 3 + E 3 ) 2Bh 

-(£3 + ^3) 2£ 2 

2Ei 2E 2 B 3 + E 3 J 



(134) 



To determine the covariant coordinates x^ = x^ — 0^ A v (|1 1 1|) , we have to fix a gauge. A 
natural gauge choice in the present context is the "static gauge" BqA^ = 0, so thatQ 



E=-diA {x), B = VxA(x). 



(135) 



The metric can be brought into diagonal form using the diffeomorphism x^' = x^ + ^(x) 
with £f(x) = 2p6(0,0,0,A (x)), which gives 



h'..„ = + d^ u + d u ^ = p0 



flU 



-2U(x)t--p 



( -(S 3 + E 3 ) 



V 
/ \ 
000 
000 

\ h(x) J 





<B 3 + E 3 ) 





\ 






[B 3 + E 3 ) 

B 3 - 3E 3 J 



(136) 



where we used (|125|) to write E 3 = B 3 + ■h h, hence 



B = E 








(137) 



2 alternatively one can also impose e.g. Aq = 0, but then the Ai become x°-dependent 
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For h(x) = 0, this has precisely the form of the metric (|121 ) with Newtonian potential 

2U = p9 (B 3 + E 3 ) = -p(26 d 3 A + ^h) (138) 

including leading relativistic corrections, while for h ^ it agrees with the Newtonian limit 
but the last term in (II 361) violates the relativistic corrections. 



To determine explicitly for given U(x), we act with 83 on ( |138|) and combine the 
result with the Bianci identity for B 

= V-B = V-E - ^d 3 h(x) = -AA (x) - ^d 3 h(x) . (139) 



This gives 

d 3 ^P- = -2d 2 3 A A.y,) =(-<)r ; -i <)-;-• ^u, (no) 

pti 

which can be solved for A as 



A (x) = y e f d 'y G ( x - y^ u (y) • ( 141 ) 

Here G(x — y) is a 3-dimensional propagator, (— d 3 + d\ + d 2 )G(x — y) = S^^x — y). The 
(static) Bianci identity for E 

= Vx£ = Vx J B + —(d 2 h(x), -d x h{x), 0) (142) 

then determines the conserved current 

J = VxB = -^-(d 2 h, -dih, 0), V-J = 0, (143) 
AO 

so that B = V x A can be solved for A. Therefore for an arbitrary given potential U(x), 
we can indeed find A^ corresponding to a static metric fluctuation h^ u which reproduces 
the Newtonian potential U(x). There is some freedom in the solution of A Q (|140|) , and h(x) 
is (almost) determined by ( |139| ). Note that even though a preferred direction x 3 is singled 
out through 9^ v and x°, this merely amounts to preferred coordinates x^ for the desired 
geometry. 

Let us consider the vacuum case AU(x) = in more detail. By integrating ( P-38|) , we 
can obtain a Ricci-flat solution with A^ = 0, h(x) = 0, 

Aq(x) = \ [ dsU(x\x 2 ,s) + H(x\x 2 ) (144) 
& Jo 

which solves AA = provided (df + dl)H = — | -£-gU(x)\ x a =0 . 
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Now consider the case with non-vanishing mass distribution AU(x) = 47rGm(x) ^ 
in a region of space near the origin. Then the presence of a (hyperbolic!) propagator in 
( |141|) implies that Aq will not be harmonic even in regions where the mass density m(x) 
vanishes. This in turn implies e.g. through ( |139| ) that h(x) 7^ 0, and the leading relativistic 
corrections to Newtonian gravity are not correctly reproduced. In other words, while it is 
possible to reproduce e.g. the Newtonian potential U(r) ~ - for a point mass, it implies in 
the electromagnetic picture a nontrivial charge density which is not localized at the origin, 



leading to a violation of Ricci flatness. This is in accord with the results of Section 3.1.1 



We conclude that the consideration of nontrivially embedded branes in matrix models 
with extra dimensions is required in order to obtain a gravity theory which reproduces the 
leading relativistic corrections of General Relativity in the presence of masses. This is in 



accord with the results of Section |2.6| that N = A SUSY is required at the quantum level, 
leading to the D = 10 IKKT model and hence to embedded branes. Since the embedding 
degrees of freedom can be viewed as scalar fields, their quantization is straightforward, and 
expected to be well-behaved. 



4 Solution with spherically symmetric Poisson struc- 
ture 

In this section, we discuss an exact but unphysical solution of the tree- level e.o.m. (|50D , 
in order to illustrate nontrivial geometries and the covariant formulation. The solution 
is unphysical, because the induced gravity action is not taken into account; this will be 
explored elsewhere. We start from the covariant e.o.m. (^T|) for the Poisson structure 9^ v 

G rn (x) V 7 0- x = -G^{x) 6-Jd 7 (j + e- 2<T G^ d 7 rj{x) (145) 

and look for a solution 9~l which is static and spherically symmetric. This Ansatz is 
actually not appropriate in order to look e.g. for a Schwarzschild-like solution; for that 
purpose one should presumably look for a deformation of the flat Moyal-Weyl solution, 
where breaks rotational invariance. Nevertheless, finding a nontrivial exact solution of 
( |145|) is certainly instructive. 

To illustrate the case of nontrivially embedded branes, consider a 4-dimensional brane 
A4 A C M 5 in the matrix model (|23|), with Cartesian coordinates x a = (x^, (p) given by the 
semi-classical limit of the matrices X a . We also use the radial variable r 2 = x\ + x\ + x\ 
and the Euclidean time r = £4. This leads to the following spherically symmetric closed 
2-form 

6T 1 = J 2) + /(r, r)dr A dr (146) 

where oj^ = sin(0)d9 A d(j) can be interpreted as field of a magnetic monopole on S 2 , which 
is singular at the origin. This will define a spherically symmetric metric G fJLV {x) if the 
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induced metric g^ v {x) on Ai A is spherically symmetric. Hence we consider an embedding 
function </> = <p(r), so thatf^ 

g^x) = Spy + d^(f)d u (f) 

ds 2 g = r 2 (d9 2 + sin 2 (9)d<p 2 ) + (1 + (f>' 2 (r))dr 2 + dr 2 (147) 

or 

g rr = 1 + <p' 2 (r), g TT = l, gee = r 2 , g w = r 2 sin 2 (6») . (148) 
Since 9~l = sin^), this gives 

ds 2 G = r- 2 (d9 2 + sm 2 (9)d<p 2 ) + f(r, rf ((1 + (j)' 2 (r)y l dr 2 + dr 2 ) . (149) 

The effective metric G^ u = e u G^ u is 



G /0,r) \v/l + (/)' 2 (r) 



where 



and 



l^| 1/2 __ /(r,r) 



|<u| 1/2 rV 1 + <l>' 2 (r) 



(150) 
(151) 



_ l/ l + /2 (r) A 
4 G ^"2l /V) +r J' (152) 



We could now define 



f{r,r) = r 2 , 
y/l + ,2 (r) = (1- 

j2/i 2 



R 2 (l-^) = r 2 (153) 

which reproduces the (Euclidean) Schwarzschild metric. However, we are not free to choose 
the embedding function (p(r), which must satisfy an e.o.m. which for the bare matrix model 
is given by = (fW|), modified by the quantum effective action, cf. (|104|) . Therefore 
the above metric is only an illustration how nontrivial geometries may be realized. We 
leave this for future work, and proceed to give an illustrative solution only for the case of 
flat embedding = resp. D = 4. 



13 



A seemingly more general 0(r, r) could be reduced to the above through a redefinition of r 
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Flat embedding <fi = 0, or D = 4. Now consider the purely 4D case with <f> = 0. The 
effective metric (|150|) then becomes 

G ab = — ^— (dd 2 + sin 2 (d) dip 2 ) + r 2 f(r, r) {dr 2 + dr 2 ) (154) 
f(r,T) 



where 

Computing the Christoffel symbols for the metric ( |154|) gives 

It = r- 1 + i/(r,r)- 1 aj(r,r) = -f; r = f; T 
f; r = ^/(r,r)- 1 a T /(r,r) = -f; r = -f;, 

f r = -G rr fd r f-\ f T = -G rr fd T f- 1 . (156) 
The covariant Maxwell equations (|145[) for the r component is then 

G rr - f r / + G Tr T r TT f - G rr f T rT f = e- 2a G TT 9 Tr d r rj(x) - G rr 6^d r a (157) 
which gives 

f(r,r)- 1 d r f(r,r) = -2f 2 r 3 (158) 
using (|155|) . Similarly, the r component of (|145 ) 

-G TT d T f + f T f-G rr f T rr f + G TT f r rr f = -e- 2 ° Grrf^r)- 1 d rV (x) + G TT f(r,r)d T a (159) 

gives 

f{r, r)- l d r f{r, r) = -e" 2<T r 4 d r r/(a;) + 5 r a = (160) 
which implies / = f(r). Together with ( [L58| ) we obtain 

/(r) = - ? 4=. (161) 
V r 4 + c 

To make contact with the standard notation, denote 



i? 2 = /-i = v ^4+^ 9 



2- ^ /. C 



Then R 3 dR = r 3 dr, hence 



~dR~^ 



(163) 



and we obtain 
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Therefore the effective metric ( 154) becomes 



ds% = R 2 {d9 2 + sin 2 (6)dcp 2 ) + J I - dr 2 + -^rdR 2 (165) 

which is flat for c = 0. This metric, in particular the radial dependence is quite strange. 
However this is not surprising, because it is a solution of the "bare" equations of motion 
only, without taking into account the induced gravitational action. Therefore this serves 
merely as an illustration how nontrivial solutions can arise. Since the (9~J we used is far from 
the Moyal-Weyl plane, the results of Section || do not apply, and there is no contradiction 
with the fact that ( |165| ) is not Ricci-flat. Indeed, note that 

^' = ~ h (166) 

which is far from the Moyal-Weyl case where e~ a = const. In particular, even the solution 
f = r~ 2 with flat is very different from the Moyal-Weyl plane. This shows how the 
same geometry may be realized in different ways. These different realizations however will 
be distinguished once the induced gravitational action is taken into account, which includes 
in particular an action for the dilaton-like field a [|TT| . 



Cartesian coordinates. To clarify the above solution, reconsider the spherically sym- 
metric symplectic form ( |146[ ) for = 0. The most general rotationally invariant antisym- 
metric tensor in 3+1 (Euclidean) dimensions has the form 

6m = x i /( r > r )> = e ijkXk g(r, t) (167) 

which is actually also invariant under 5*0(3). The corresponding 2-form 

6~ l = f(r, t) Xidrdxi + g(r, r)eijkXkdxidxj (168) 

is closed if and only if 

g (r) = r~ 3 , (169) 

recovering (|146|) 



6 1 = f(r, t) Xidrdxi + r 3 SijkXkdxidxj (170) 
for an aritrary function /(r, r). The corresponding effective metric is 

Goo = e£s»e£ = f*f 

Gu = e^s^ + e^e^^r-H^ + xMf-r- 6 ). (m) 

For / = r -3 , we obtain 



G -1 

Z 1 " r 4 





' 





8u t 



(172) 



33 



so that G^ v reproduces the flat solution found above ( |165|) for c = 0. Note again that the 
corresponding (9~J is very different from the Moyal-Weyl case where = const. 

We conclude that in order to obtain realistic metrics such as the Schwarzschild- metric, 
different nontrivial embeddings must be used, solving the equations of motion derived from 
the combined bare action plus induced gravitational action. 



5 Symmetries and conservation laws 



The basic matrix model (23) is invariant under the D-dimensional Poincare group, consist- 
ing of translations 

X a -> X a + c a , c a eR (173) 
and rotations resp. Lorentz transformations 

X a ^A a b X b , A b G SO(D — 1, 1). (174) 

These symmetries lead to conservation laws according to Noethers theorem, which are 
elaborated below for the case of translations; see also for a related discussion. Adapting 
a standard trick, we consider the following non-constant infinitesimal transformation X a — > 
X a + 5X a for 

5X a = {X\[X a ,6 b ']}g w (175) 

where e b is an arbitrary matrix, and g a b = 5 a b or g a b = rjaj,. As elaborated in Appendix C, 
this leads to 

^SS YM = -Tre c \X\ f a ' c '}g aa ,g cc , (176) 

for arbitrary e a , where 
fab = [X \X c ][X\X c '] gccl + [X b ,X c }[X a ,X c '}g cc , - ±g ab [X d ,X c }[X d ',X c '}g dd ,g cc , (177) 



is the matrix - "energy-momentum tensor" . Since ( pL76| ) vanishes on-shell, the conservation 
law 

[X a ,f a ' c ]g aa , = (178) 

follows. This can of course also be checked directly using [X a , [X b ,X a '}}g aa , = 0. Moreover, 
since it is a consequence of a symmetry of the action, this will survive quantization in the 



form of a Ward identity. Indeed it is easy to check that (|175|) defines a measure-preserving 



vector field on the space of matrices X a , so that ( |178| ) also holds under the matrix path 



integral i.e. upon quantization; there will be additional terms in the presence of matter 
or in correlators. Note that the indices of the "tensor" T ab range from 1 to D, including 
transversal components. A covariant form of these conservation laws and their physical 
meaning in the context of gravity remains to be elaborated. 
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A very similar conserved energy-momentum tensor was obtained previously in p4|.|55 
in the context of NC gauge theory on the Moyal-Weyl quantum plane. In that case, it 
was possible to find a suitable gauge invariant version of T ab which satisfies a standard 
conservation law ||54|| . The present result is somewhat different since (|178|) is obtained for 
NC spaces with general 6^(x), involving also components which are transversal to the 
brane. Moreover, the meaning of gauge invariance versus locality is somewhat different 
(and not entirely clear) in the present context; for example, U(l) gauge transformations 
are now interpreted as symplectomorphisms. In any case, a similar "local" version of 
( p.77| ) involving Wilson lines might help to clarify its interpretation. An analogous energy- 



momentum tensor in the context of the BFSS matrix model was also found in [56 



6 Discussion and outlook 

We present in this paper a general framework for studying emergent gravity in the context 
of Yang-Mills type matrix models, on generic noncommutative branes embedded in R D . 
The basic message is that the dynamics of fields on the brane is governed by an effective 
metric in the semi-classical limit, which depends both on the embedding and the Poisson or 
noncommutative structure on the brane. The resulting geometry is dynamical, governed by 
the matrix model and its induced effective action which includes in particular the Einstein- 
Hilbert term. Therefore Yang-Mills matrix models contain some type of gravity theory. 
The results of [||] are thus generalized to a richer class of geometries, setting the stage for a 
systematic exploration of the physical properties of the models. This necessity to consider 
nontrivially embedded branes in higher dimensions is shown by a detailed analysis of the 
Newtonian limit of the D = 4 model, which does not correctly reproduce the relativistic 
corrections to the Newtonian limit. 

Matrix models such as the IKKT model therefore provide a simple and transparent 
mechanism for gravity, which arises from fluctuations of the basic matrix degrees of free- 
dom, along with nonabelian gauge fields. While the IKKT model was proposed originally 
as non-perturbative description of IIB string theory flT2| , [57| 1, the progress in this and related 



works shifts the emphasis towards the consideration of general noncommutative branes and 
geometries, which promise to provide the physically relevant backgrounds. They appear to 
be simpler and more natural in this context than classical spaces and geometries, the essen- 
tial difference being the effective metric which involves the noncommutative resp. Poisson 
structure. Similar considerations should apply also to time-dependent matrix models such 
as the BFSS model |2l|. 



There are some important differences to General Relativity. The essential point is that 
the metric is not a fundamental degree of freedom, but arises effectively as described above. 
This leads to important simplifications for the quantization: first, the issue of gauge fixing 
is much simpler, involving degrees of freedom which can be viewed as scalar and gauge 
fields in a NC background. Second, it is not the Einstein-Hilbert action which is quantized, 
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rather the matrix model action, which is similar to a Yang-Mills action. This allows to 
compute e.g. the one-loop effective action in a straightforward way, which boils down to 
computations in a NC gauge theory or the use of standard heat-kernel expansions under 
certain conditions. Most remarkably, in the case of maximal supersymmetry (i.e. the IKKT 
model in D = 10) the model can be expected to be finite, leading to the identification of 
the Planck scale with the scale of N = 4 SUSY breaking. This suggests that the IKKT 
model may provide a well-defined quantum theory of fundamental interactions including 
gravity. 

Remarkably, emergent NC gravity appears to provide a mechanism for avoiding the 
cosmological constant problem, which is explained in the case of 4-dimensional branes. 
Again, the full significance can only be judged once near-realistic solutions are found and 
understood. Here the compactification of higher-dimensional NC branes as indicated may 
turn out to be important, which is motivated also from particle physics, providing a mecha- 
nism for gauge symmetry breaking and fermionic zero modes. A full discussion of emergent 
gravity in such cases is a challenging subject for future work. 

This paper contains only semi-classical considerations. These are the leading terms 
in a systematic expansion in O^, which should be elaborated eventually. This can be 
achieved using the Seiberg-Witten map P7| , which allows to systematically re- write a 
noncommutative (gauge) theory in terms of a commutative one. While this was used in [|J 
to obtain the semi-classical limit of the nonabelian gauge fields in emergent gravity, it is not 
part of the definition of the model: it is simply - by definition - a natural way to extract 
the physical content of a NC model. In principle, the quantization should be done on the 
level of the matrix model, and its effective action can then be interpreted in a commutative 
language. For example, the issue of UV/IR mixing is resolved here not through the Seiberg- 
Witten map but through its proper interpretation in terms of gravity @,|TT[. At least in 
the case of (softly broken) N = 4 SUSY, one may hope to resolve similarly the issues of 
unitarity and Wick rotation. All this clearly requires much more work. 

Let us summarize the main arguments supporting emergent NC gravity as described by 
.D- dimensional matrix models: 

• The models do describe some gravity theory on 4-dimensional NC branes, since matter 
couples to a universal metric (up to possibly density factors). Gauge fields and gravity 
are naturally unified. 

• The class of geometries is rather rich in the case of models with D > 4. 

• The geometry is dynamical, governed by an effective action which includes the 
Einstein-Hilbert term at the quantum level. The quantization is likely to be well- 
defined, at least for the IKKT model. 

• Flat space is a solution even at the quantum level, without fine-tuning 

• The models are extremely simple, without any classical-geometric prerequisites. 
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This certainly describes a very promising theory of gravity, the main missing item being 
the analog of the Schwarzschild solution. This requires to consider nontrivial embedding 
as shown here, and is complicated by the fact that the quantum effective action is required 
at least at one loop. 
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Appendix A: Some identities 



The following is an important identity for Poisson tensors: 



-d^9^' u - 26 v ' v p- l d vl p 



(179) 



noting that 2p 1 d u p = d v 9 m> 9 h, hence 



3>r) = o. 



180) 



On-shell vanishing of T^: For our restricted class of metrics, the above identity ( |180| ) 
together with \G pu \ l l 2 = pe CT fl54]) implies 

r> = -\G pa \- l l 2 d u (&»\G pu \ 1 ' 2 ) = --e-°d v (G^p) 
1 



—<r°d v { P 9 w e^g^(x)) 

-e- a 9 w 'd v {9^'g^{x)) c '= 







;i8i) 



using the e.o.m. ([|8|) for X M . This can also be seen from fl44|). Therefore the equations of 
motion for are equivalent to = 0. From the point of view of General Relativity, this 
would be interpreted rather as a gauge-fixing condition. This is not the case here due to 
the constrained class of metrics. 
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Appendix B: Derivation of the covariant e.o.m. 



Consider 



where 

Using (|36|) 
where 

we can write 



G ab Tl 



,d p {CPo^G ab ). 



\fG ah 

CT = e~° G"" = \x)e vv ' \x) g^(x) 
9n'p'{x) = e~ u g^^(x), 

= -G 1V G p5 6~l {d 1 G pu + d v G pl — d p G 
— ( d^G pu + d y Gp^ — d p G^ u 



G -yv G pS 9 -i drr G pi/ 



71/ 



e^d^Gpv 



= -G^(d,G pS )9-tG pv 

= -G lv 9~l Gp,(d^gp 5 ')e 55 ' + e^gps'd^ 55 ') 

= -G^G pv g m d^ - G^Gp^gp, + G™^" 1 

since 9 1P := G cr, G ps 8~l is antisymmetric. The Jacobi identity gives 

G^~g m d^ = e^'-g^rg^'d^ 

= 0^~g lW ~g m (6 p W,6^ + 6 P W^' P ) 

= 6^~g jW 6^~g m d 7 9^ + (6^6^~g m yg Yr)l d^'f 

= e^g^'g.pd^'' - G™g m d,6 pp 



hence 



G^~g m d^ 



-e^~g^'~g m d^' 



Finally, observe that 
G^d^-G^d^ 



G 1V d 7 (d^d v <j>) - G ir} dp(d j( pd v <p) 
dp4> G^ d^d v cf) + G 1V dpd^d^ 

dp^G^d.d^-^dpSg^ 



(If 



12) 

13) 

14) 
15) 



(186) 



(187) 



(If 



dpd v (j) 
(189) 
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hence 

d^g m = l -G^ d p 5g yv + G 1V d 7 d v( j> (190) 

and 

& m d l9m = ~G™ d^g 7V + d P 4> G in d,d ri 4> . (191) 

Therefore 

G 1V d 7 g m = \& m dtfTr, + e~° d p <P G 7V <W - G^g m d,a + l -G^~g in d p o (192) 
and we obtain 

G^~g m + &W»d 7 ~g m = U^-g^'-g^d^' + V G™d y g m 

+0™ (e-° G^d 7 d v 4> - G^~g m d,a + d p a) 
= p » (drf(x) + e- CT d p cf> G^d^d v (j> - & m ~g m d^a + 2f,(x) d p a) 

using the scalar function 

fj(x) = \&^g pu = -^g^g^ (193) 

which satisfies 

d 7 fj(x) = U^'-g^'-g^ + ^GTOyg^ . (194) 
Putting all this together, we obtain 
G^CaOV^- 1 = G^d^l-G^T^Q-l-^Q-l 

= G pv (G^~g m d^ + G^d^ - f%l 

= G pv (e^id^ix) + e~ CT d p <P G^d 7 d^ - &mg m d,a + 2fj(x) d^afj - f%, 

This holds identically, i.e. it characterizes the constraint of the metric. 

Now we take into account the equations of motion 1^ = ( |181| ) and = G^ u d p d v 4> = 
(f46|) , which hold in the special coordinates x p defined by the dynamical matrices (this is 
why the above is non-covariant). We can then rewrite this as 

G™(x)%6£ = G p J^ (d^(x)-G^g m d J a + 2f}(x)d^ 

= G pv 6^ (d^(x) + 2fj(x) d p a) - G^Q~l d y a . (195) 

Using n = e 2a fj(x) = \9^' g p i u ,6 uv ' g pu (§9|), we obtain 

G^{x) V 7 (e ff B~l) = e a G pu e pl (d 7 fj{x) + 2fj{x) d 7 a) = e~ a G pu 6^ d 7 r)(x) (196) 

which is flT5|) resp. ([5(]). This is the covariant form of the equation of motion, independent 
of the choice of coordinates. 
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Appendix C: Derivation of ([ITS]) 



We use here a short-hand notation where double upper indices are understood to be con- 
tracted with 5 a b or r]^. Then 

= —SSyu = Tr[5X a ,X b ][X a ,X b ] 

= Tr[{X c , [X a ,e c }},X b }[X a ,X b ] 

= Tr[X c [X a ,e c ] + [X a , e c ]X c , X b ][X a , X b ] 

= Tr(^X c [[X a ,e c },X b }[X a ,X b } + [[X a , e c ], X b ]X c [X a , X b ] 
+ [X a ,e c ][X c ,X b ][X a ,X b ] + [X c ,X b ][X a ,e c }{X a ,X b ]^ 
= Tr ({X c , [[X a , e c ] , X b ] } [X a , X b ] + [X a , e c ] { [X c , X b ] , [X a , X b ]}^j 
= Tr ( - {X c , [X b , [X a , e c }}}[X a , X b ] + [X a , e c ]{[X c , X% [X a , X b ]}^j 
= Tr(±{X c , [e c , [X b ,X a ]]}[X a ,X b ] + [X a ,e c ]{[X c ,X% [X a ,X b ]}) (197) 
for arbitrary e a . Using Tr({A, [B, C]}C) = Tr([A, B]C 2 ) this can be written as 
= Tr^-^[X c ,e c ][X a ,X b }[X a ,X b } + [X a ,e c }{[X c ,X b },[X a ,X b ]}^ 

= Tr^e c [X c , [X\ X b ] [X a , X b ]] - e c [X a , { [X c , X b ] , [X a , X b ] }]) 

= Tre c [X a ,f ac } . (198) 
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